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ABSTRACT 

We compare the spatial distribution of stars which form in hydrodynamical simulations 
to the spatial distribution of the gas, using the Q-parameter. The Q-parameter enables 
a self-consistent comparison between the stars and gas because it uses a pixelated 
image of the gas as a distribution of points, in the same way that the stars (sink 
particles in the simulations) are a distribution of points. We find that, whereas the 
stars have a substructured, or hierarchical spatial distribution {Q ^ 0.4 — 0.7), the 
gas is dominated by a smooth, concentrated component and typically has Q 0.9. 
We also find no statistical difference between the structure of the gas in simulations 
that form with feedback, and those that form without, despite these two processes 
producing visually different distributions. These results suggest that the link between 
the spatial distributions of gas, and the stars which form from them, is non-trivial. 

Key words: stars: formation - ISM: structure - star clusters: general - methods: 
numerical 


1 INTRODUCTION 

The physics of star formation results in stars grouped to¬ 
gether in regions which exceed the mean density of the 
Galactic disc by several orders of magnitude ( BlaaraJ 196J; 


iLada &: Ladall20o3 : iPorras et al.ll2003 : iBressert et al.ll2010ll . 

Furthermore, there is growing evidence that star formation 


is hierarchical from the interstellar medium fISM) down 

to sub-PC scales (iHovld 1 19531: IScalc 

1 19851: lEfremov 

1 19951: 

lElmegreen et al.ll2006l: Bastian et al. 

2007l:lKruiissenl 

2 OI 2 I), 


and that there is no preferred spatial scale for a given star 
formation ‘event’. 

Understanding this complex hierarchical picture of star 
formation requires analysis of the spatial distribution of gas 
from the ISM down to the giant molecular clouds (GMCs) 
from which stars form, to the substructure of the clouds and 
then of the spatial distribution of the stars themselves. In 
particular, do stars exhibit the same spatial distrib ution as 
the ISM (e.g. as argued by iGouliermis et all l2014l . for the 
NGG346 star-forming region) and if so, does this provide 
information on the physical processes from the ISM down 
to individual star formation? Answering this question re¬ 
quires the study of star-forming regions at the earliest pos¬ 
sible stages, since the gas structure is likely to be disrupted 
by feedback and the stellar structure by dynamical inter¬ 
actions. Either or both of these processes can in principle 
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erase any spatial correlation between gas and stars on short 
timescales. 

In order to address these questions, analysis of 
both gas in the ISM and GMGs, and stars in young 
regions, must be undertaken using a self-consistent 
method. In recent years, structural analysis of stars 
in young regions has b een performed using the pow¬ 
erful Q-parameter fe.g. ICartwright fc WhitworthI l2004l: 
BastiaB_et_alJ_ 200^ Schrnmaetal 20081 : ^nchez fc Alfarol 


J3agtiajy^t_aJjj2UU^ljbcnmeia^t^i.l uuua : isancnez fe Aitarc 
2 OO 9 I : Gregorio-Hetem et al.ll20lJ b which combines informa¬ 


tion on the minimum spanning tree of a distribution with 
the typical separation between the points in the distribu¬ 
tion. This technique has also been developed to study the 
gas distribution in images, by appropriately weighting the 
flux from pixels in an image to create a d istribution of points 
ibomax. Whitworth fc Gartwrightll201ll '). 

It is not clear to what extent and for how long the spa¬ 
tial distribution of stars should follow the same distribution 
as the gas from which they form, and this has yet to be 
addressed observationally. However, detailed hydrodynami¬ 
cal simulations of star formation provide both information 
on the spatial distribution of gas, and on the distribution 
of stars (sink particles). In this paper, we exam ine the hy¬ 
drody namical simulations of star formation by IPale et aU 
kOlJl and measure the spatial distribution of both the gas 
(from pixelated images) and the distribution of the sink par- 
ticles, usin g the Q-parameter. The suite of simulations by 
IPale et ^ include control runs without feedback from pho- 
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toionisation and stellar winds, and so we also search for dif¬ 
ferences in the gas distributions for clouds that are influ¬ 
enced by feedback mechanisms, and those that are not. 

The paper is organised as follows. In Section [2] we de- 
cribe the implementation of the Q-parameter on gas distri¬ 
butions, in Section[3]we present our results, in Section|4]we 
provide a discussion and we conclude in Section [5] 


2 METHOD 

The simulations of lPale et al.l use a Smoothed-Particle Hy¬ 
drodynamics (SPH) code to model the evolution of GMCs 
of a range of sizes and masses and are seeded with Burgers 
turbulence such that their initial virial ratios are either 0.7 
or 2.3. The formation of stars is followed using sink par¬ 
ticles llBate et al.lll995ll which, in the simulations analysed 
here, have accretion radii of 0.005 pc and can be regarded 
app roximately as individual stars. 

IPale et ^ allow their simulations to evolve until three 
objects exceed 20 M© at which point each calculation is 
forked into a control run and a feedback run. The control 
runs evolve purely hydrodynamically as before, while the 
feedback runs are impacted by the ionising radiation and/or 
stellar winds of the massive s tars, model l ed res p ective ly by 
the algorithms presen ted in IPale et ^ (l2007l . l2012l l and 
IPale fc Bonnelll (l2008l ). All runs are continued for as near 
as possible to 3 Myr to evaluate the effects of pre-supernova 
feedback on the clouds and clusters. 

For the analysis presented here, pixelated column- 
density maps are constructed by drawing a pixel grid over 
the simulation, placing each SPH particle on the grid and 
column-integrating through its smoothing kernel onto all 
grid cells whose centres it overlaps. Since the maps presented 
here have a lower resolution than the local resolution of the 
hydrodynamic simulations, particles which are too small to 
overlap the centre of any grid cells have their mass smeared 
out over the area of the cell in w hich they lie. _ 

We use the Q-p arameter (ICartwright fc WhitworthI 
I 2 OO 4 I : [Cartwright|l2009l l to quantify the spatial distribution 
of both stars and g as in the hydrodynamical simulations 
of iPale et ^ (120141 . The Q-parameter is determined by 
constructing a minimum spanning tree (MST) of all of the 
points in a distribution and then dividing the mean MST 
branch length, fti by the mean separation between points, s: 



( 1 ) 


Petermining Q for the stars (i.e. sink particles) in the 
simulations is tri vial, but determining Q for the ga s 
is more involved JCartw right. Whitworth fc Nutt^I 2 OO 6 II . 


Ihomax. Whitworth fc CartwrighT i 20111 1 provide a method 

for converting the flux distribution in a pixelated image 
into a distribution of points, from which the Q-parameter 
can then be calcul ated. We refer the interested reader to 
I Lomax et alJ (|201ll l for full details of the method, including 
examples of its use on synthetic images, but b riefly sum¬ 
marise the method here. In the simulations of IPale et ^ 
(|2014h . we use gas column density as the flux. 

For an A^pix = X x J array of pixels, the total flux 



Figure 1. Dependance of the Q-parameter on the number of par¬ 
ticles in a distribution. From top to middle, the lines indicate a 
smooth, concentrated distribution with power law density profile 
n oc (solid purple), (dashed raspberry), (dot- 

dashed pale blue), (dotted orange), r® (solid magenta); and 

from bottom to middle the lines indicate a substructured fractal 
distribution with fractal dimension D = 1.6 (solid black), 2.0 
(dashed red), 2.6 (dot-dashed green), 3.0 (dotted dark blue). The 
error bars indicate the interquartile range from 100 identical re¬ 
alisations of the same cluster for a given number of points. 


received from the pixels, Ftot is: 

( 2 ) 

In order to convert the flux d istribution into a distribution 
of points, iLomax et al. then define the flux quantum 

as 

(3) 

We start by choosing a pixel, TZij at random. If Fij ^ AT, 
then the flux in that pixel is reduced 


Fij —>■ Fij — AT 


( 4 ) 


and we place a point at nj + Armd, where Vij is the centre of 
pixel TZij and Armd is a small random displacement an order 
of magnitude smaller than the pixel size, to prevent the final 
point distribution from having a gridded appearance. 

If Fg < AF, then we consider a patch of n pixels (see 
fig. 1 in lLomax et al n is increased until the flux from 

the patch is equal to, or exceeds AF: 


Fti —patch ^ Fi, ^ AF (5) 

n —patch 


The flux from each pixel is then uniformly reduced 


Fij —>■ Fij 



AF \ 

^71 — patch / 


( 6 ) 


and a point is placed at position rpnt which is equal to the 
weighted centre of the removed flux, plus a small random 
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displacement: 


rpnt = X] 

n —patch 


Fr,.- 


n —patch 




( 7 ) 


This process is repeated A^pix times, whereafter the total flux 
is now zero. 

We now have a distribution of A^pix points with 
which we calculate the Q-parameter for the gas distri¬ 
bution. The Q-parameter provides a measure of the de¬ 
gree to which a distribution is s ubstr u ctured, or concen - 
trated. Cartwright fc WhitworthI (l2004l l. rCartwrightI (l2009l l 
and iLomax et al.l ( 201 il l provide calibration data for syn¬ 
thetic models, so for a given Q-parameter one can assign a 
fractal dimension (in the substructured case), or the radial 
density profile exponent (in the smooth, concentrat ed case). 

However, as mentioned bv iLomax et al.l ll201lf) . the Q- 
parameter varies slightly depending on the number of points 
in the distribution. This is usually not a problem when com¬ 
paring e.g. the outcome of simulations to real star-forming 
regions, as the numbers of stars are similar (between 100 
- 1000). However, when constructing a point distribution 
from pixels, in order to properly sample the gas distribu¬ 
tion we find at least 64 x 64 = 4096 pixels are required, but 
128 X 128 = 16384 pixels are desirable. 

In Fig. [T]we show the dependence of Q on the number 
of points, N. The coloured lines from top to middle indi¬ 
cate a smooth, concentrated distribution with power law 
density profile n oc ® (solid purple), ® (dashed rasp¬ 


berry). 


(dot-dashed pale blue), r (dotted orange). 


(solid magenta). From bottom to middle, the lines indi¬ 
cate a substructured fractal distribution with fractal dimen¬ 
sion D = 1.6 (solid black), 2.0 (dashed red), 2.6 (dot-dashed 
green), 3.0 (dotted dark blue). The error bars indicate the 
interquartile range from 100 statistically identical realisa¬ 
tions of the same cluster for a given number of points. Q 
tends to more extreme values with increasing Ai if a distri¬ 
bution is highly substructured, or highly concentrated, but 
remains roughly constant for more moderate distributions. 
In order to check that these results are not simply due to 
sampling too few points, we created each dataset with 10“^ 
points and then randomly chose N from this distribution 
and found very similar results. 


3 RESULTS 


We an alyse five sets of SPH simulations from IPale et ^ 
dioli). These simulations follow the evolution of giant 
molecular clouds as they form stars without feedback 
from photoionisation and stellar winds (‘control runs’) 
and simulations in which feedback is switched on (‘dual 
feedback runs’). In earlier work dParker fc Paid 120131 : 
IParker. Pale fc Ercolanol l2015h . we showed that the pres¬ 
ence (or not) of feedback can effect the long term spatial 
evolution of a star-forming region and here we investigate 
whether the spatial distribution of the gas is influenced by 
feedback, and if the spatial distribution of the gas follows 
that of the stars. 

The results from the control Run I simulation 2.2 Myr 
after the epoch when the calculation is split into parallel 
control and dual-feedback incarnations are shown in Fig. (2] 


gas. We also show the positions of the sink particles (the 
white points), which have Q = 0.72, suggesting a slightly 
substructured distribution. In Fig. 2(b) we show the distri- 


In Fig. 2(a) we show the surface density distribution of the 


button of points from the pixel distrib ution of column den - 
sity, determined using the method from iLomax et ahl ll201lll . 
This gas distribution has Q = 1.01, suggesting a smooth, 
concentrated distribution. 

We show the results at the s ame timeste p of th e dual¬ 
feedback Run I calculation from IPale et ^ ll2014ll which 
includes the effects of photoionisation feedback and stel¬ 
lar winds in Fig. |3] The sink particles (shown by the white 
points in Fig. 3(a) I have Q = 0.49, suggesting a substruc¬ 
tured distribution. From inspection of Fig. 3(b)[ the dis¬ 
tribution of gas also appears substructured; however, the 
Q-parameter is 0.88, indicating a smooth, concentrated dis¬ 
tribution. 

The results for all five sets of simulations are sum¬ 
marised in Table [T] (all five sets have very similar morpholo¬ 
gies). We present the Q-parameter for a pixelated map of 
the gas distribution for each simulation using 64 x 64 pixels, 
and 128 x 128 pixels (differences between the two are min¬ 
imal). In the 128 x 128 pixel case, we also present the Q- 
parameter for cold gas only (i.e. ionised particles have been 
removed), and for cold gas above a column density threshold 
of > 1 X 10“^gcm“^ (this process sets the column density 
of around 25per cent of the pixels to zero). Finally, we give 
the Q-parameter for the sink particles, and the number of 
sink particles in each simulation. 

For one simulation set (Run I) we also present the Q- 
parameters for the simulations at the point before feedback 
is switched on. The gas distribution in this simulation is 
(as we might expect) similar to the end of the control run 
(Qgaa = 1.02), but the Q-parameter for the sink particles 
is similar to the value at the end of the dual-feedback run 
Qsinks = 0.42. This is likely due to the simulated cluster 
being substructured early on, and the control run erasing 
some of tha t substructure due t o dynamical evolution (as 
discussed in IParker fc Dalel[201^ . 

Two main results are apparent in our analysis. Firstly, 
the Q-parameter for the sink particles is systematically lower 
than the Q-parameter of the gas distribution. [Note that 
this is not due to the differences in the numbers of parti¬ 
cles used to determine the Q-parameter, as highlighted in 
Fig- □ this effect is only important if both Q-values are 
in the same regime (i.e. substructured, or concentrated)]. 
Typically, Qgas ~ 0.9 (indicating a smooth, concentrated 
distribution), whereas the sink particles tend to have Q- 
parameters between 0.4 - 0.7 (indicating a substructured, 
or hierarchical distribution). 

Secondly, and perhaps most strikingly, differences be¬ 
tween the gas distributions in the simulations with and with¬ 
out feedback are minimal in most cases (despite the appar¬ 
ently substructured gas distribution for the runs in which 
feedback is switched on). If we compare the value of Qgas 
for the control run simulations, the range of values is 0.86 - 
1.01 for the 128 x 128 pixel maps. The dual feedback simu¬ 
lations range from Qcoid gas,i 28 x 128 = 0.78 — 0.92. However, 
individual simulations often exhibit almost identical values 
for Qgas (e.g. the feedback runs from simulations J, UF and 
UP have Q values higher than, but very similar to, the con¬ 
trol runs). Furthermore, differences in Q of a s little as 0.1 
should not be taken as being significant (e.g. IParker et al.1 
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(a) Raw column density image of the gas and sinks (white dots) in the (b) Particle point representation of the gas only in the control 
control Run I simulation. For the sink particles, Q = 0.72. Run I simulation at a resolution of 128 x 128 pixels; Q = 1.01. 

Figure 2. The spatial distribution of stars and gas in the control Run I simulation from iDale et alj <l2014h 2.2 Myr after the point where 
the calculation is forked into separate dual-feedback and control runs. In panel (a) we show the surface density distribution the gas, and 
the sink particles are shown by the white points. In panel (b) we show the pixel point distribution for the gas. The sink particles have 
Q = 0.72 and the gas distribution has Q = 1.01. 



(a) Raw column density image of the gas and sinks (white dots) in the (b) Particle point representation of the gas only in the dual- 
dual-feedback Run I simulation. For the sink particles, Q = 0.49. feedback Run I simulation at a resolution of 128 x 128 pixels; 

Q = 0.88. 

Figure 3. The spatial distribution of stars and gas in the dual-feedback Run I from lPale et al.l ll2014l) . 2.2 Myr after the enabling of 
photoionisation and winds from the 0-type stars. In panel (a) we show the surface density distribution of the gas, and the sink particles 
are shown by the white points. In panel (b) we show the pixel point distribution for the gas. The sink particles have Q = 0.49 and the 
gas distribution has Q = 0.88. 
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Table 1. Measured Q-parameters for the gas distributions in the simulations from lPale et alj ||2014|V The columns display the run ID, 
whether feedback was switched on, the Q-parameter for the gas distribution using 64 x 64 pixels, the Q-parameter for the gas distribution 
using 128 X 128 pixels, the Q-parameter using 128 x 128 pixels for cold gas only, the Q-parameter using 128 x 128 pixels for cold gas 
with a column density above 1 X 10“^gcm~^, the Q-parameter for the sink particles, and the number of sinks. 


Run ID 

Feedback 

Qgas,64 X 64 

Sgas,128 X 128 

Qcold gas, 128 x 128 

Qcold gas, flux limited,128 x 128 

Qsinks 

-^sinks 

I 

Off 

1.03 

1.01 

1.01 

0.96 

0.72 

186 

I 

Dual 

0.87 

0.88 

0.83 

0.81 

0.49 

132 

I 

Before 

1.02 

0.99 

0.99 

0.95 

0.42 

44 

J 

Off 

0.91 

0.89 

0.89 

0.89 

0.49 

578 

J 

Dual 

0.87 

0.90 

0.91 

0.92 

0.70 

564 

UF 

Off 

0.88 

0.86 

0.86 

0.84 

0.77 

66 

UF 

Dual 

0.83 

0.85 

0.88 

0.89 

0.49 

93 

UP 

Off 

0.92 

0.88 

0.88 

0.86 

0.49 

340 

UP 

Dual 

0.87 

0.88 

0.92 

0.88 

0.64 

343 

UQ 

Off 

0.91 

0.86 

0.86 

0.88 

0.70 

48 

UQ 

Dual 

1.07 

0.83 

0.78 

0.77 

0.45 

77 


1 


0.8 


0.6 

o 

JO 

E 

0.4 


0.2 

0 

0 0.2 0.4 0.6 0.8 1 

sbar 

Figure 4. The rfi — s plot llCartwrightlbOO^ for different mor¬ 
phologies containing N = 10000 points. We create 100 realisations 
of each morphology. Anti-clockwise from the bottom: fractals with 
fractal dimension D = 1.6 (the black crossed circles), D = 2.0 (the 
red hashtags), D = 2.6 (the green dotted circles), D = 3.0 (the 
dark blue crosses); then radially smooth clusters with power law 
density profiles n oc (the magenta plus symbols), n oc r~^'^ (or¬ 
ange asterisks), n oc r~^'^ (blue plus symbols), n oc ® (rasp¬ 
berry hexagrams) and n oc ® (purple compressed squares). 
We overplot the rh and s values for the pixel-point distributions 
for the gas in the simulations. Control runs are shown by the open 
symbols, dual feedback runs are shown by the filled symbols. Fi¬ 
nally, values for rh and s corresponding to a Q-parameter of 0.8 
are shown by the solid black line. 



l2014h and based on the apparent visual structure in the feed¬ 
back runs, we would expect Q values of less than 0.8. Given 
the apparent differences between the column-density images 
in F igures [2] and [3l t his re sult seems counterint uitive. 

ICartwr^htl |20^) and iLomax et ahl ll201ll) provide a 
further diagnostic of the underlying spatial distribution in 


relation to the Q-parameter by plotting fh against s. In 
Fig. |4]we show the expected values of fh and s for distribu¬ 
tions of = 10 000 points with various morphologies, sam¬ 
pling 100 versions of each morphology. Anti-clockwise from 
the bottom: fractals with fractal dimension D = 1.6 (the 
black crossed circles), D — 2.0 (the red hashtags), D = 2.6 
(the green dotted circles), D — 3.0 (the dark blue crosses); 
then radially smooth clusters with power law density pro¬ 
files n oc. (the magenta plus symbols), n oc (or¬ 
ange asterisks), n oc (blue plus symbols), n oc ® 

(raspberry hexagrams) and n oc (purple compressed 

squares). Q = 0.8 is shown by the solid line. 

We plot the values of fh and s for the pixel-point dis¬ 
tributions of the gas in each of our simulations; the control 
runs are shown by the open symbols, and the dual feedback 
runs are shown by the filled symbols. Run I is shown by the 
triangles. Run J by the circles. Run UF by the stars, UP by 
the squares and Run UQ by the diamonds. The values of fh 
and s for the control runs sit mainly around (but not on) 
the parameter space of the smooth, concentrated profiles 
with n oc (the blue plus symbols). These fake pro¬ 

files typically have Q-parameters between 0.9 - 1.0, which 
is roughly consistent with the values in the simulations (0.86 
- 1.01). The runs with feedback have more scatter, and most 
sit roughly between the smooth, concentrated profiles with 
n oc and n oc ° (orange asterisks). The n oc ° 
profiles have Q-parameters between 0.8 - 0.9, again similar 
to the simulation values of 0.78 - 0.91. 


Despite Q, fh and s indicating a smooth, concentrated 
distribution, the pixel-point distributions for the gas in the 
simulations with feedback appear (at least to the eye) sub¬ 
structured. In order to test whether this is due to a failing 
of the Q-parameter in determining the structure of complex 
distributions, we perform two simple Monte Carlo tests in 


an attempt to mimic the spatial distribution in Fig. 3(b) 


First, we create a smooth, concentrated distribution of 
10 000 points with radius 10 pc and a radial density profile 
n oc as shown in Fig. 5(a) This has a Q-parameter 

of 1.1. We then remove all points within 2 pc of the origin, 
and move them to the outskirts of the distribution (thereby 
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Figure 5. The effects on the spatial distribution of a smooth, concentrated cluster with density profile n oc ® when removing centrally 
located points. In panel (a) we show the original distribution, which has Q = 1.1. In panel (b) we have removed everything within 2 pc 
of the centre and placed it at larger radii, following the underlying radial profile, and Q = 0.75. In panel (c) we have removed everything 
within 5 pc and placed them at larger radii, and Q = 0.60. 


maintaining the same number of points and the same den¬ 
sity profile) and Q — 0.75 (Fig. 5(b) I. We then repeat this 
process, but remove everything within 5 pc of the centre, 
and Q — 0.60 (Fig. 5(c) I. This is intended to mimic in an 
approximate fashion the sweeping up of gas and clearing 
out of bubbles, which is the main visible effect of feedback. 
However, we keep the distribution of points smooth. Whilst 
the result is very similar in appearance to the gas distribu¬ 
tion from the dual-feedback run (Fig. 3(b) I, it clearly has 
a very different underlying spatial distribution according to 
the Q-parameter. 


Secondly, we create a ‘broken ring’ of 2000 points, which 
has Q — 0.3. This is designed to represent more faithfully 
the broken and irregular inner walls of the feedback-driven 
bubbles. We then embed this in a uniform field of a further 
1000 points (with Q — 0.7), representing the smooth and 
largely undisturbed background gas, and then finally, place 
a further 3000 points in a smooth, concentrated sphere (with 
an n oc ® density prohle - Q = 1.7). This material rep¬ 
resents the very dense clumps of gas found in the feedback 
runs. This final distribution, which incorporates all three 
components, is shown in Fig. [6] and has an overall Q = 0.9. 
Clearly, the concentrated clump in the centre of the distri¬ 
bution is dominating the overall Q-parameter. 


In order to test this in the SPH simulation data, we 
take the pixel-point distribution for the dual-feedback run 
shown in Fig. 3(b) and remove the region of highest density 
centred on {3,4} pc. This distribution is shown in Fig. [T] 
When we remove this region, the Q-parameter is reduced 
from 0.88 (indicating a smooth concentrated distribution) 
to 0.76, which is in the slightly substructured/uniform field 
regime. The remaining ‘background’ of points, which are 
tracing low-density gas in the simulation, are still contribut¬ 
ing to the Q value, hence the distribution still does not ap¬ 
pear very substructured. 


Our interpretation is therefore that the gas distribu¬ 
tion in the dual-feedback runs - whilst visually appearing 
to be very different from the control runs - is actually very 



x/pc 

Figure 6. Three different spatial distributions overlaid within the 
same field. We show a broken ring of 2000 points within a uniform 
distribution of a further 1000 points. In the centre is a centrally 
concentrated sphere of a further 3000 points with an n oc ® 
density profile. In isolation, the broken ring has Q = 0.3, the 
uniform field has Q = 0.7 and the sphere has Q = 1.7. The 
combined distributions have an overall Q-parameter of 0.9. 


similar, and both are dominated by a smooth, concentrated 
component. 


4 DISCUSSION 

We set out to examine the extent to which the structure of 
the gas and stellar distributions in a set of hydrodynamic 
simulations are similar using an objective quantitative cri¬ 
terion: the Q-parameter. This is in contrast to many other 
studies where the surface density of either the stellar mass 
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Figure 7. The gas pixel-point distribution for the dual-feedback 
Run I simulation 2.2 Myr after feedback is enabled (Fig. |3(b)| l, but 
with the pixel-points from the most dense region of gas removed. 
Removing this region of points reduces the Q-parameter from 0.88 
to 0.76. 


or the star formation rate are compared with the surface 
density of the gas, which can only offer qualitative conclu¬ 
sions. 

If our results can be taken at face value, they suggest 
that the spatial distribution of stars is very different to that 
of the gas from which they form. In one sense, thi s is not sur- 
*risin g; the initial conditions of the simulations of lPale et al.l 
2014h are smooth, concentrated gas clouds, whereas stellar 
distributions, particularly in the early stages of star forma¬ 
tion, are expected to be highly substructured. However, once 
the cloud has evolved and formed stars in a substructured 
distribution, we might expect the gas from which they form 
to also exhibit substructure. This is not the case, and the gas 
distribution remains dominated by a smooth concentrated 
component. This is particularly evident in Run I where the 
largest subcluster is fed by a set of smooth accretion flows 
which do not themselves fragment into clumps or stars. 

The presenc e (o r not) of feedback was shown in 
IPale et al] (l2012l 'l and IParker fc PaI3 (l201,3l ') to have a min¬ 
imal effect on the spatial distribution of stars at the end of 
the feedback calculations, wh en compared with the control 
runs - all th e simulations fromIPale et al.l l|2012h . [Pale et ahl 
(l2013ll and IPale et ^ (l20l4l form stars with a substruc¬ 
tured dist ribution, and other authors report similar conclu- 
sions (e.g. lBatdl20i^ : lOirichidis et al.l[2012l : iKruiissen et al.1 
I 2 OI 2 I I. From visual inspection, it appears that feedback re¬ 
sults in a very different spatial distribution for the gas (com¬ 
pare Figs. 2(b) and 3(b) I. However, most of the gas remains 
in a dense, concentrated distribution(s) in these pixel-point 
maps. The difference in fact is in the relative positions of the 
stars and gas. In the control simulations, the densest con¬ 
centrations of gas and stars are spatially coincident. In the 
feedback simulations, the accretion flows feeding the clusters 
are destroyed or deflected and the clusters and their envi¬ 
rons are cleared of gas. Much of the gas is swept up into 


relatively dense concentrations, but these form stars ineffi¬ 
ciently, so that the main effect of feedback is to spatially 
separate some of the stars from dense gas (this can readily 
be seen by comparing the positions of the sink particles to 
the gas in Fig. 3(a) I. However, in most of the runs, even in 
the presence of feedback many of the sink particles are still 
coincident with the dense gas in 2P projection, even if they 
not still genuinely embedded in three dimensions. This is 
particularly true of Run J, where large quantities of dense 
gas remain projected at the location of most of the stars, 
even though in three dimensions the main concentration of 
stars has been largely cleared of gas. This effect blurs out 
the dynamical influence of feedback when seen in projection. 

One caveat is that the Q-parameter may not be an ideal 
method for measuring the spatial distribution of gas, as we 
first need to conver t a pixelated image into a distribution 
of points. However, iLoniax et al.1 (l201ll l show that the un¬ 
derlying spatial distribution is always recovered in synthetic 
datasets, and in the Monte Carlo experiment shown in Fig. [6] 
we also recover the underlying spatial distribution. 


5 CONCLUSIONS 


We measure the spatial distributions of gas and stars in 
hydrodynamic al simulations of star formation using the 
Q-pa rameter dCartwright fc Whitworthll2004l : iLomax et al.1 
I 2 OI j l. We compare Q for simulations which form under the 
influence of feedback from photoionisation and stellar winds, 
and corresponding control run simulations. Our conclusions 
are as follows. 

(i) The spatial distribution of the gas is different to 
that of the stars in all simulations. The gas has a spatially 
smooth, concentrated distribution (Q ~ 0.9), whereas the 
stars have a substructured distribution (Q ~ 0.4 — 0.7). 

(ii) The presence of feedback clears out cavities, or bub¬ 
bles in the simulations, which are visible to the eye. However, 
statistically these distributions are similar to the control run 
simulations without feedback. The reason for this is that in 
all simulations, a significant dense, concentrated gas compo¬ 
nent dominates the spatial distribution. 

(iii) The combination of points (i) and (ii) suggests that 
a direct link between the spatial distribution of young stars, 
and the gas from which the stars form should not necessarily 
be expected in observations. 

These results also highlight the pitfalls in trying to mea¬ 
sure quantitatively something previously judged by eye. The 


column density images shown in Figures 2(a) and 3(a) and 


the corresponding pixel-point distributions could scarcely 
look more different, yet the Q-parameter analysis reports 
that they are statistically indistinguishable. The root cause 
of the similarity is that they are both dominated by the same 
kinds of structures, the visual difference being due to how 
those structures are arranged. 

There is growing evidence in the literature that sug- 


scale tefremov 

199f 

lElmeereen et al.1 l200fil: iBastian et al.1 

I 2 OO 7 I: IKruiissen 

2012 

h. Because of this, it is tempting to link 


different regimes together, such as the spatial distribution 
of gas in the ISM with the spatial distribution of stars in 
clust ers (e.g. iGouliermis et al.l l20l4 iGresorio-Hetem et al.1 
I2OI5II . Our results suggest that - whilst star formation in 
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numerical simulations also produces a hierarchical spatial 
distribution of the stars - the distributions of gas and stars 
may be linked in a highly non-trivial way. 
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